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ABSTRACT 



BuckminsterfuUerene (Ceo ) has recently been detected through its infrared emission bands in the interstellar medium (ISM), includ- 
ing in the proximity of massive stars, where physical conditions could favor the formation of the cationic form, CJ^. In addition, 
CggWas proposed as the carrier of two dilfuse interstellar bands in the near-IR, although a firm identification still awaits gas-phase 
spectroscopic data. We examined in detail the Spitzer IRS spectra of the NGC 7023 reflection nebula, at a position close (7.5") to the 
illuminating B star HD 200775, and found four previously unreported bands at 6.4, 7.1, 8.2, and 10.5 yum, in addition to the classical 
bands attributed to polycylic aromatic hydrocarbons (PAHs) and neutral Cgo . These 4 bands are observed only in this region of the 
nebula, while Cgo emission is still present slightly farther away from the star, and PAH emission even farther away. Based on this 
observation, on theoretical calculations we perform, and on laboratory studies, we attribute these bands to Cgg. The detection of Cg^ 
confirms the idea that large carbon molecules exist in the gas phase in these environments. In addition, the relative variation in the 
Ceo and CJ„ band intensities constitutes a potentially powerful probe of the physical conditions in highly UV-irradiated regions. 



1. Introduction 

The mid-infrared (mid-IR) spectrum of galactic and extragalac- 
tic objects exhibits band emission (strongest at 3.3, 6.2, 7.7, 
8.6, and 11.2 fim) attributed to carbonaceous macromolecules, 
i.e., polycyc lic aromatic hydrocarb ons (PAHs, see recent state 
of the art in iJoblin & Tielensll20TTh . In addition to PAH bands, 
IR emission b ands at 7.0, 8.5, 17.4, and 19.0 um ha ve recently 
been reported jCami et al l20ialSellgren et alJ2010h . and found 
to match qui te closely th e IR active bands of buckminster- 
fuUerene (Ceo, iKroto et aD[l985i) . a cage-like carbon molecule. 
Carbonaceous macromolecules, including PAHs, carbon clus- 
ters, or fullerenes, are believed to play a fundamental role in 
the physics and chemistry of the interstellar medium (ISM), and 
their infrared signatures are commonly used as a tracer of physi- 
cal conditions. Nevertheless, Ceo is the only molecule belonging 
to this family, which has been specifi cally identified in the ISM. 
In the NGC 7023 reflection nebula, ISellgren et al] ( |20IO|) have 
shown that Ceo is predominantly found in the regions closest 
to the star. In that part of the nebula, UV irradiation is high 
(above 10"* times the interst ellar standard radiation field), an d 



PAH molecules are ioniz ed dRapacioli et alj|2p05 
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^^^^^ ■ ■20I2h. i f not destroyed (iBerne & Tielens 

20I2t iMontillau d et al.ll20I2l) . One could therefore expect Cgg 



to be present in these regions. 

iFoing & Ehrenfreundl d 1 9941) reported evidence of interstel- 
lar CgQ based on the detection of two diffuse interstellar bands 
(DIBs) at 9577 and 9632 A, however this identification is 
still questioned considering that no spectrum of Cg,, could be 
recorded yet, in conditions appropriate for DIB identification, 
i.e., in gas phase and at low temperature. The IR s pectrum of 
CgQ was measured in a rare-gas matrix by iFulara et al. Cl993i) 
and wa s found to exhibit two bands at 7. 1 and 7.5 ^um. lKern et alj 
(l20I2h have performed new spectroscopic measurements and 



suggests that the latter band is du e to Cg^, whereas the authors 
attribute a band at 6.4 yum to C^q. iMoutou et"al1 (1 1 9991) derived 
upper limits on the abundance of Cg,, in NGC 7023 based on 
not detecting the 7.1 (and possibly wrongly attributed 7.5 /im) 
bands. So far, there has been no observational evidence of any 
7. 1 or 6.4 jjm bands in astronomical sources. 

Looking carefully at the Spitzer data of NGC 7023, we found 
four emission bands, at 6.4, 7.1, 8.2, and 10.5 fim, which are only 
present in the regions closest to the star This also corresponds 
to a region where Ceo emission is strong. A natural carrier to 
explain these bands is Cg^, and this assertion is supported by 
spectroscopic arguments that we discuss hereafter. 



2. Observations 

NGC 7023 was observed with the short-low (SL) and long-low 
(LL) modules of the I nfraRed Spectrog raph (IRS. iHouck et al.l 
l2004l onboard Spitzer (i Werner et al.ll2004<) in spectral mapping 
mode. The spectral resolution of IRS is A/AA = 60 - 130. The 
slit width (comparable to angular resolution) is 3.6" for the SL 
and 10.5" for the LL module s. Data reduction was performed 
with the CUBISM software (ISmith et al.l l2007h and consisted 
in cube assembling, calibration, flux correction for extended 
sources and bad pixel removal. From the LL data cube we ex- 
tracted maps integrated in the H2 line at 17.0 fim, the PAH band 
at 16.4 yum and the Ceo band at 19.0 ^m. The maps ar e presented 
in Fig. 1 in a fashion similar to ISellgren et al] (l20I0l) . From the 
SL data cube, we extracted two significant spectra corresponding 
to different regions in the nebula (Fig. 1). The extraction regions 
were 3x5 pixels (1.6"/pixel) to improve the signal-to-noise ratio. 
Position 1 corresponds to the cavity of atomic gas between the 
star and the molecular cloud at an angular distance of 2 1" (~ 0.04 
pc) from the star HD 200775 (q'=21:01:34.8, ^=+68:10:07.3). 
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Fig. 1. Left: False-color image of the NGC 7023 nebula such as the one presented in lSellgren et al.l (1201 Ol) . obtained from integrating 
different components in the Spitzer-IRS LL spectral cube. Red is the emission integrated in the Ceo 19yum band. Green is the 
emission of the PAH 16.4 fim band. Blue is the emission integrated in the H2 (0-0) S(0) 17.0 jum band. The white rectangles indicate 
the regions where the IRS-SL spectra shown in the right panel hav e been extracted. The black dashed rectangle depicts the IRS-SL 
field of view. The white dashed rectangle shows the region where [Sellgren et al] ( 1201 Ol) extracted their spectrum. Right: Spectra at 
positions 1 and 2 in the image. The spectrum at position 1 has been shifted up and scaled down to allow easier comparison to the 
spectrum at position 2. Error bars are not shown here but are comparable to the width of the line (see Fig. 2). The red lines indicate 
the four newly detected bands attributed to C^q. The DFT calculated spectrum (with scaled wavelength positions) is shown as a bar 
graph in blue. 



Position 2 (ff=21:01:37. 3, 5=4-68:09:55.3) is the closest position 
to the star in the IRS spectral cube (7.5", ~ 0.015 pc from HD 
200775, see Fig. 1). These two positions correspond to regions 
where Cgo emission is strong in the 19.0 jum band (Fig. 1). For 
both positions the signal-to-noise ratio is very good and ranges 
between 100 and 200. 

3. Observational results 

The spectra (Fig. 1) have in common that they are dominated by 
bands at 6.2, 7.7, 8.6, and 11.2 lu m, which are a ttributed to vi- 
brational modes of PAH molecules (lTielensl2005l) . In addition to 
these bands, the spectrum at position 2 shows several bands that 
are absent in the rest of the nebula. These are at ~ 6.4, 6.6, 7.0, 

7.1, 8.2, and 10.5 jum (Fig.[T]). These bands are seen in several 
pixels of the IRS cube and well above the instrumental errors 
(Fig. 2). To derive the precise parameters for these bands we fit 
them using Gaussian profiles (Fig. 1) and splines to subtract the 
underlying emission due to the wings of PAH bands. The po- 
sitions, widths, and intensities of the bands are given in Table 
1. Th e 7.0 fim band has been attributed to Cgo ( ISellgren et alj 
I2OIOI) . The 6.6 yu m band has recently been attr ibuted to (pos- 
sible) planar C^a dGarcia-Hernandez et aDl201 ll) . The 6.4, 7.1, 

8.2, and 10.5 fim bands have not been observed or discussed 
yet. Since HD 200775 is a B star, only low ionization potential 
atoms should emit in fine structure lines (e.g. [CII], [SIII] etc.), 
and these species do not have lines in this spectral range, so we 
exclude contamination by fine structure lines. These four new 
bands seem spatially correlated, i.e. all of them only appear in 
the regions close st to the star, which suggests a common car- 
rier. As shown bv lSellgren et al.l (l2010l) . Ceo is also found only 
close to the star. Still, the four new bands only appear in the 
regions that are closest to the star, while Fig. 1 demonstrates that 
Ceo emission is more extended. This suggests that the four new 



bands are carried by a species that is a product of the photopro- 
cessing of Ceo , an obvious carrier being Cg^. In the following 
section we provide the spectroscopic arguments that support this 
observational evidence. 



4. Spectroscopy of C^g 

The only IR spectro scopic data of has been obtained in 
rare gas matrices by iFulara et"an { 1993h and more recently by 
iKern et aTl (|2012|) . Two bands at 7.1 and 6.4jum seem defini- 
tively attributed to C^g, based on these experimental studies. 

Theory is another approach to obtaining spectroscopic data, 
and density functional theory (DFT) in particular has been 
shown to be effective and accurate for calculations on neu- 
tral C eo (IChase et al ] 1 19921: lFabiarilll996t llglesias-Groth et al.1 
l20Tlh . However, an additional theoretical problem is that upon 
ionization Ceo is known to undergo Jahn-T eller (JT) distortion 
dChancev & O'BrienlfT997tlBersukerll2006h . Neutral Ceo has //, 
symmetry, and is a closed-shell system, and its ground elec- 
tronic state is totally symmetric and nondegenerate. Upon ion- 
ization, the ground electronic state of Cg,,, in the adiabatic ap- 
proximation, is five-fold degenerate and has h,, symmetry. These 
electronic states, degenerate at the symmetric configuration, split 
when the symmetry is broken, and the Jahn-Teller theorem pre- 
dicts that there must be lower energy extrema, some of which 
must be minima, at distorted geometries. The closeness of these 
electronic states means that the adiabatic approximation may 
break down in this situation. Nevertheless, standard harmonic 
vibrational spectra can be computed in the adiabatic approxima- 
tion around the new minima, and they will still be approximately 
coiTect if IR-active modes have a negligible component along 
the directions involved in the Jahn-Teller effect. This is the case 
for Cg,), as discussed in Appendix lAl where the reader can also 
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Table 1. Properties of fitted bands for position 2 



Position FWHM 



Intensity 



I(P W m-' sr- 



Species 



Assignment 



6.24 


0.13 


1.02 


PAH 


6.43 


0.09 


3.68 


^60 


6.61 


0.13 


3.41 


C24 


7.02 


0.07 


2.92 


Ceo 


7.13 


0.11 


2.82 




8.10 


0.10 


0.91 


^60 


8.49 


0.14 


1.83 


Ceo 


8.62 


0.28 


0.85 


PAH 


10.53 


0.11 


0.22 


^60 



Table 2. Computed IR-active bands of D5,, C 
modes with an intensity above 0.1 km.mol ' 



gQ. Only IR active 
are shown. 



Symmetry 


Unsealed 


Scaled 


Intensity 


degeneracy 


frequency 


wavelength* 






(cm-') 


ifim) 


(km.mor') 


El (2) 


1607 


6.40 


94.6 


Azd) 


1599 


6.43 


5.5 


Azd) 


1488 


6.91 


5.6 


El (2) 


1482 


6.94 


2.4 


El (2) 


1434 


7.17 


123.8 


El (2) 


1353 


7.60 


33.3 


El (2) 


1333 


7.71 


1.4 


A2(l) 


1250 


8.23 


20.0 


El (2) 


1249 


8.23 


39.5 


El (2) 


1196 


8.60 


1.4 


Ajd) 


989 


10.4 


6.6 


El (2) 


981 


10.5 


22.6 


El (2) 


795 


12.9 


2.0 


El (2) 


778 


13.2 


6.8 


El (2) 


763 


13.5 


0.5 


A2(l) 


592 


17.4 


3.5 


El (2) 


588 


17.5 


0.3 


A2(l) 


553 


18.6 


19.8 


El (2) 


412 


25.0 


12.1 


El (2) 


364 


28.2 


5.0 


A2(l) 


361 


28.5 


0.6 



using an empirical scaling factory = 0.9725 

find a more detailed description of the Jahn -Teller effect and 
how it affects our calculations. 

We performed our DFT calculations using the hybrid B3LYP 
exchange-correlation functional and the 4-3 Ig Gaussian basis 
set. This combination is known to yield reasonably accurate vi- 
brational frequencies, after scaling with an empirical factor x 
to account for the overestimation of the frequencies. The rela- 
tively small basis set limits the accuracy of band intensities , but 
based on the thorough tests performed bv lGalue et al.l (1201 ih . we 
expect that neglecting dynamical JT effects is the leading limi- 
tation on vibrational band accuracy even at this level of theory. 
The distorted geometry of minimum energy was obtained by op- 
timizing with no symmetry constraints, and turned out to have 
D^y symmetry. This is only very slightly distorted. Such a small 
change in geometry implies a correspondingly small change in 
interatomic forces, and it explains the qualitative similarity of the 
computed spectra of and Ceo- Of course, many more modes 
are IR-active in distorted due to the lowering of symmetry 
and correspondingly relaxed selection rules. The properties of 
the most intense IR active bands of Ci, are given in Table 2. 




8.2 8.4 8.6 8.8 9 
Waveknglh dim) 



0.4 10.6 
Waveknglh dim) 



Fig. 2. Fits (dashed line) to the observed bands (diamonds) in 
position 2. The continuum has been subtracted. Red solid lines 
denote emission attributed to C^, blue to Ceo , green to C24 , and 
black to PAHs. 



5.1. 



C+o in NGC 7023 

Identification of C+,/n NGC 7023 



As mentioned in the previous section, it is necessary to correct 
the calculated frequencies by an empirical scaling factor. In the 
case of PAHs, this is usually done by comparing the computed 
frequencies to the ones measure d in the laboratory at low tem- 
perature in rare gas matrices (see lBauschlicher & LanghoftI 19971 
for a case-study on PAHs). Therefore we use the IR absorption 
spectrum of Cg^ that was measured in Ne matrices by Kern et al. 
(2012) to calibrate our DF T calculations. Th e two bands defini- 
tively attributed to Cg^ by dKern et al.ll20I2l) are found at 1550 
and 1406 cm ' and correspond to the strongest bands, predicted 
by theory at 1607 and 1434 cm"' . This implies a respective scal- 
ing factor of 0.9645 and 0.9805. We adopt an average value 
X - 0.9725. The corrected positions are given in Table 2 and 
the resulting spectrum is shown in Fig. 1. After scaling, the five 
strongest IR bands fall at wavelengths of 6.40, 7.17, 7.60, 8.23 
and 10.50 jjm. Four of these are very close (within 2%) to the 
positions of the four new bands detected in NGC 7023 (6.43, 
7.13, 8.10, 10.53 fim). The observed match is very good consid- 
ering that other factors are expected to affect the band positions, 
in particular ba nd shift due to anh armonic coupling in hot emit- 
ting molecules (iJoblin et al.lll995h . The nondetection of the 7.6 
/im band in the observations is not surprising, since it is most 
likely hindered by the strong PAH emission at the same position 
(Fig. 1). Based on the observational (Sect. O and spectroscopic 
arguments presented in this Letter, we therefore argue that there 
is strong evidence for the presence of C^q in NGC 7023. Further 
discussion of the relative band intensities would require an emis- 
sion model. Attempts to build such a model for Cm have been re- 
ported by Sellgren et al. (2010) and then by Bernard-Salas et al] 
il2012h . However, the authors could not account for the observed 
relative band intensities. This is likely due to uncertainties in the 
intrinsic IR band strengths that were obtained by DFT calcula- 
tions. The same certainly holds for Cg,, for the reasons given in 
Sect.|4] Still, the detectio n of Cg^ strongly sup ports the idea that 
Ceo is in the gas phase dSellgren et al. 201 Oh a nd not in solid 
phase as initially suggested by Cami et alT ( 2010h . 
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5.2. Abundance 

Using a simple energetic consideration, we can derive an esti- 
mation of the Cggabundance. We write that C(,o and C^q re- 
lax in the IR all the energy they absorb in the UV, which is 
true if IR emission is the only relaxation channel. We can fur- 
ther assume that the UV absorption cross-section of neutral and 
cationic species are simil ar, which is the case for PAH s at en- 
ergies higher than 10 eV dCecchi-Pestellini et aDl2008h . In this 
case, the ionization fraction of Ceo can be directly derived from 
the integrated IR flux in the Ceo and Cg,)bands. Using the values 
in Table 1 we find 4.7 x 10"^ Wm^^sr"' for Ceo and 7.7 x 10"^ 
Wm^^sr"' for C^,,. However, these values have to be corrected by 
the emission at longer wavelengths. In the case of Ceo , we there- 
fore add the emission in the 17.4 and 18.9 yum bands, measured 
in the IRS2 LL spectrum for position 2 (respectively, 2.6 and 
5.2 xlO"^ Wm^^sr^'). This leads to a total value of 12.5 x 10"^ 
Wm^^sr ' for the total IR flux. In the case of C^q, we do not 
consider emission at longer wavelengths in line with DFT cal- 
culations, which predict much weaker bands. Thus, we obtain 
that the ionization fraction of Ceo is 38%. If we consider a max- 
imum abundance of Ceo in NGC7023 of 1.7 x 10"'' of the el- 
emental carbon (iBerne & Tielensl[2012h . this implies an abun- 
dance of ~ 1.0 X 10"* of the elemental carbon abundance for 
C gQ. This is a facto r at least 10 lower than the value derived by 
iFoing & Ehrenfreun d ( 1994) for the diffuse interstellar medium, 
assuming that Cg^is the carrier of the two DIBs at 9577 and 9632 

A. Since the DIB identification is still not firmly established, it 
is diflicult to discuss this discrepancy further. 

5.3. Taking the next step, C^q as a tracer of physical 
conditions 

ISellgren et al.l (1201 Oh report the detection of the Ceo band at 7.0 
fim in NGC 7023. In their spectrum, the bands of C^q are not 
present. Since this spectrum was taken in a region farther away 
from the star (Fig. 1) this suggests that the Ceo emission ex- 
tends farther away from the star than C^,, emission. These rela- 
tive variations can be attributed to the photochemical evolution 
of Ceo resulting from the competition between ionization by 
UV photon and recombination with electrons. Assuming one can 
quantify these two processes, and also have a proper description 
of the photophysics in these systems, this would allow using the 
ratio between the IR flux of Ceo and as a tracer of local 
physical conditi ons. There already ex ists such an approach us- 
ing PAH bands dGalliano et alj|2008l) . but only empirical laws 
can be used since the PAH population is poorly characterized. 
Furthermore, information on the Ceo and bands is unique 
for regions of high UV radiation fields in which other molecular 
tracers may not survive. 



6. Conclusion 

After studying the mid-IR spectra of the NGC 7023 nebula, we 
have found spectral signatures at 6.4, 7.1, 8.1, and 10.5 yum, 
which we attribute to the cationic form of Ceo (Cgg). This is the 
largest cation know in space so far. has been proposed as a 
DIB carrier, and our identification supports this proposal. This is 
also clear evidence for the presence of large carbon molecules 
in the gas phase in the ISM. The detection of Cg,, in emission 
also opens the possibility of using the ratio between the Ceo and 
CgQ IR bands as a tracer of physical conditions in interstellar 
and circumstellar environments. This may prove useful in the 



framework of forthcoming infrared missions, such as JWST and 
SPICA. 
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Appendix A: Details on Cg^spectroscopy 

Neutral Ceo has //, symmetry, its highest occupied molecular or- 
bital (HOMO) is five-fold degenerate, it has h,, symmetry and is 
fully occupied, resulting is a closed-shell ground state that is to- 
tally symmetric and nondegenerate. Upon ionization, the hole in 
the h„ HOMO yields a five-fold degenerate /j,, overall electronic 
state. This undergoes spontaneous symmetry-breaking due to 
the Ja hn-Teller (JT) effect (Chancey & O'Brien 1997; Bersuke] 
I2006h . The degeneracy of the electronic state is lifted by distort- 
ing the molecule to a lower symmetry along some of its normal 
modes, which are determined by symmetry and called JT-active. 
For CgQ , JT-active modes are those of Hg, Gg, and Ag symme- 
try. The Ag modes only shift the total energy, without reducing 
symmetry. The Hg, Gg modes instead break the //, symmetry, and 
produce a multisheet adiabatic potential energy surface, with a 
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conical intersection in the symmetric geometry and extrema in 
configurations of D5,, and Dj,^ lowered symmetry. DFT predicts 
that the Dsy geometries should be the ab solute minim a, with the 
Djr ones being shallow transition states (ISaitoll2002h . However, 
high-resolution photo-electron spectroscopy seem s to hint that 
the D 31, geometry could be the real minimum jCanton et alj 
I2002h . When JT-distorted minima are deep with respect to vi- 
brational energy {static JT-effect), the adiabatic approximation 
holds in its vicinity, and standard harmonic vibrational analy- 
sis is applicable. Conversely, if equivalent minima are separated 
by negligible potential barriers, the molecule can tunnel among 
equivalent minima, mixing the near-degenerate electronic states 
(dynamic JT-effect), and the adiabatic approximation is not ap- 
plicable. The resulting vibronic states recover the full initial 
symmetry of the problem. In this case, a much more complex 
calculation, dropping the adiabatic approximation, is needed for 
accurate results. A comparable, but much simpler case of dy- 
namical JT effect is the cation of corannulene, C2oH|q, which 
can be regarded as a fragment of Ceo wi th periphera l bond s 
saturated by H atoms. This was studied bv iGalue et al.l ( 1201 Ih . 
who compared experimental infrared, multiphoton dissociation 
(IRMPD) spectra with a plain DFT harmonic vibrational analy- 
sis at the distorted geometry of minimum energy. This was ex- 
pected to be the worst possible comparison, since in IRMPD ex- 
periments, vibrational energy is increased until the dissociation 
threshold is reached. This corresponds to energy values that are 
much higher than all barriers among equivalent minima, thereby 
maximizing dynamical JT effects. Despite this, the experimen- 
tal and theoretical spectra do qualitatively agree, allowing for 
accurate band identification. Bands that are most significantly 
mispredicted, in position and intensity (but still identifiable in 
the laboratory spectrum), are those whose normal modes imply 
displacement along JT-active modes, i. e. those that move the 
molecule from a minimum in the direction of another one, or to 
the conical intersection. 

In the light of this, we did a similar analysis for CJg, finding 
the distorted geometry of minimum energy and computing har- 
monic vibrational spectra there, thereby neglecting dynamical JT 
effects. We performed DFT calculations using both the Gaussian 
version 03. d2 and NWChem version 6.1 codes, and obtained 
very nearly identical results. Optimization led to the D5,, geome- 
try. We also optimized the geometry of CggWith the constraint of 
//, symmetry, obtaining a JT stabilization energy of ~ 70 meV, 
which is consistent with previous calculations (Saitol l2002l) . The 
Dsv distorted geometry, when compared with the symmetric one, 
appears to be distorted almost exclusively along normal modes 
of Hg symmetry, with changes in bond lengths of a few mA and 
bond angles by less than a degree (maximum). 



